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ABSTRACT 


An algorithm that extracts aerosol optical depth and particle size index from satellite 
observations of upwelling radiance over clear ocean arcas was modified and studied. In 
order to examine the algorithm’s performance on a regional scale, a previously analyzed 
data set, retrieved bv the NOAA-9 AVEIRR sensor, was reprocessed. The area of the 
study was in the central Pacific Ocean during the RITS-88 crutse from 7 April to 5 May 
1988. The results were compared to those of the earlier study and used to investigate 
changes caused by modifications to the algorithm. One significant modification was to 
account for the absorption of column water vapor in the calculations. In addition, the 
correlation between aerosol optical depth and the two-term Henvevy-Greenstein phase 
function was investigated by normalizing measured radiance values. 

The results compared well with the earlier study and supported the modifications to 
the algorithm. Based on this study, aerosol optical depth was shown to be correlated to 
the Henyey-Greenstein phase function. This knowledge helped in formulating improve- 
ments to the phase function. 

It was concluded that satellite imaging and processing on a regional scale is a useful 
wav to study marine aerosols, that results can be rmproved by considering the effects of 
atmospheric water vapor absorption and other modifications to the calculations, and 


that there 1s room for 1mprovement to the two-term Henvey-Greenstein phase function. 
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I INTRODUCTION 


A. THE NATURE OF THE PROBLEM AND BACKGROUND 

As the world’s population and knowledge have grown, so too has an interest in its 
climate. The ability to study our environment and anticipate possible climate changes, 
whether caused by humans or not, has far-reaching implications for future generations. 
Their very survival may depend not only on our ability to understand any changes, but 
also on What we do about them, whether it be to adapt to the changes or attempt to 
influence our climate. How well we might adapt to changes or whether we can influence 
our climate is still not known but we have begun to understand it by studying the overall 
radiation budget of the earth. 

Prior to the advent of the satellite, estimates of radiation entering and leaving the 
earth-atmosphere system Were made using radiative transfer equations {Rao et al. 1990). 
Many factors had to be considered and incorporated into the model. Radiation 1 the 
atmosphere 1s attenuated by many processes. It 1s absorbed and scattered by molecules 
of air, aerosols, and meteorological features such as clouds, precipitation and fog. Aer- 
osols are solid and liquid particles, approximately one micron or less in diaineter, sus- 
pended in air. This includes smoke, dust, haze, and some clouds. While the results of 
the radiation calculations were helpful, there are still many uncertainties. Tor example, 
in recent years there has been increasing interest in the warming of the earth’s climate 
caused by increased concentrations of the so-called greenhouse gases, notably CO, , 
which trap outgoing long wave radiation. 

It was and still is necessary to parameterize the less well understood contributions 
of such features as aerosols and clouds. These features are inherently more difficult to 
measure than CQO, because they have shorter time scales and they arise on a local or re- 
gional scale leading to higher variability. Sources of aerosols include natural, like 
volcanic aerosols and desert dust, anthropogenic, or human made, and biological, where 
gaseous metabolic byproducts are precursors to aerosols. McCormick and Ludwig 
(1967) proposed that a buildup of aerosols could have a cooling effect on the earth. 
Charlson and Pilat (1969) concluded that particle size distribution as well as scattering 
and absorption cross sections also played a role and could lead to warming in some 
cases, especially for anthropogenic aerosols. Coakley et al. (1983) studied naturally oc- 


curring tropospheric aerosols and found that they result in global surface cooling of 2-3 


°C, an amount close to the increase in temperature due to such gases as CH, and N,O, 
or a doubling of CO,. In addition, the cooling effect of aerosols can be multiplied when 
they act as cloud condensation nuclei (CCN) resulting in the formation of cloud drop- 
lets. Clouds reflect incoming solar radiation leading to a cooling effect. Opposing this 
effect is a tendency for non-nucleated anthropogenic aerosols to absorb radiation 
thereby reducing cloud reflectance which leads to warming. Twomey (1977) and 
Twomey et al. (1984) concluded that for thin and medium optically thick clouds over 
continents, an increase in anthropogenic aerosols results in cooling. Thin and medium 
optically thick clouds are most common in the atmosphere. Thus, fluctuations in aer- 
osol concentration could play a role in the global climate with an increase in aerosol 
concentration contributing to a cooling effect. 

Furthermore, Charlson et al. (1987) has expanded on hypotheses of Shaw (1983) and 
later Lovelock (1986) that a biological feedback mechanism where marine 
phytoplankton, through production of marine aerosol precursors, play a role 1n climate 
regulation. 

Figure | is a schematic diagram of how the feedback mechanism works. In the di- 
agram the rectangles are measurable quantities and the ovals are processes linking the 
rectangles. A positive sign in the oval indicates the effect a positive change in the pre- 
vious rectangle has on the following rectangle. A negative sign has a negative effect on 
the following rectangle for a positive change in the previous rectangle. 

Charlson et al. (1987) explain that the dimethylsulfide (DMS) produced by 
phytoplankton is oxidized in the atmosphere to form a non-sea-salt sulfate (NSS-SO2>). 
These NSS-SO?z-’s act as CCN in the marine atmosphere. The argument 1s that DMS 
is produced in sufficient quantity to be the primary source of CCN in the remote, un- 
polluted marine atmosphere. Since clouds form from CCN, and clouds act to reflect 
incoming solar radiation, phytoplankton might then play a role in regulating the earth's 
climate. The greatest uncertainty in the hypothesis is what effect any lowering of the 
earth’s surface temperature and reduction of solar irradiance below the clouds would 
have on the phytoplankton and DMS production. 

Once satellites became available they were quickly tasked to help provide data on 
radiation or factors that influence radiation (Rao et al. 1990). For example, marine 
aerosols are atmospheric constituents that appear to play an important role in the 
overall radiation budget, as was shown. Hindman et al. (1984) describe experiments that 
establish a positive relationship between aerosol variations and satellite detected bright- 


ness Variations in the marine boundary layer. In fact, direct measurements of marine 
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Figure 1. Biological feedback in the Earth’s radiation budget (Charlson et al. 


1987): rectangles are measured quantities, ovals are linking processes. 


Plus or minus signs indicate the effect a positive change in a previous 
rectangle has on the following rectangle. L 


aerosol characteristics could be provided by ship or aircraft. While this information is 
valuable, it 1s expensive and difficult to obtain and only a small area of the world’s 
oceans can be sampled and studied at a time. Satellites, on the other hand, can provide 
near continuous data over large areas. This capability is desirable for evaluating the 
effects of aerosols on a global scale. Meteorological satellites, however, have been de- 
signed to measure cloud albedo while sensitivities to radiances associated with aerosol 
optical depth variations are small (Durkee 1984). Thus, in order to determine aerosol 
characteristics, 1t 1s important that the satellite data be carefully processed. 

Aerosols can also play a role in military applications. Aerosols can scatter energy 
from a laser beam propagating through the atmosphere. Lasers can be used for satellite 
communications with submarines. The narrow beam of the laser 1s harder to detect than 
the normal radio communications now in use and 1s capable of higher data transmission 
rates to deeper depths. As submarines become more important as multimission plat- 
forms, it becomes more important to have secure communications with them (Painter 
1989). In addition, line-of-sight communications operating in the near infrared can uti- 
lize aerosol scatter in the marine boundary laver to provide over-the-horizon ranges of 
30-300 miles (Mooradian 1981). 

Further investigations of aerosol characteristics are clearly needed. The purpose of 
this thesis is to better utilize satellite data in order to more accurately measure acrosol 
characteristics so that their effect on the overall radiation budget of the earth and in 


other applications such as communications may be quantified. 


B. METHODOLOGY 

Durkee et al. (1991) describes an algorithm whereby multispectral radiance data 
obtained from the Advanced Very High Resolution Radiometer (AVHRR) onboard the 
U.S. National Oceanic and Atmospheric Administration (NOAA)-9 satellite is used to 
obtain estimates of aerosol optical depth and particle size index over clear ocean areas. 
By utilizing different satellite channel sensitivities, the particle size index can be meas- 
ured and the optical depth calculated. 

Modifications to the algorithm have been made. For example, Mahony (1991) 
studied and quantified the effects of column water vapor absorption on aerosol particle 
size index. This information can be used to improve the performance of the algorithm. 
In addition, a modification to account for variable earth-sun distance on solar irradiance 
was made. A change in the computer code to correct a mathematical error and an ad- 


justment to an over correction for ozone absorption were the final two modifications. 


To study the modified algorithm on a regional scale, a data set used by Benedict 
(1989) was reprocessed. The results were compared to those of Benedict’s (1989) earher 
study. It was expected that all the modifications would improve the results with the 
water vapor correction being the most significant. An important extension of this pro- 
cedure was to study the correlation between opiical depth and the two-term Henyey- 
Greenstein (HG) phase function used in the calculations. This was accomphshed by 
studying a plot of normalized measured radiance versus scattering angle for various 
particle size indices. The purpose of doing this was to make optical depth independent 
of scattering angle; these variables might be related through the phase function. If there 
is a correlation, improvements to the phase function and therefore calculated optical 
depth are indicated. 

This methodology was chosen because the Benedict (1989) data set was available 
and provided a baseline from which to measure improvements. In addition, useful 1n- 
formation about the HG phase function had not vet been extracted from the Benedict 
(1989) data set. Once improvements to the algorithm are verified and permanently in- 
corporated into the computer code, new data sets can be studied with the goal of better 


understanding the role of DMS on the earth's climate. 


C. OBJECTIVES 

The specific objectives of this thesis are to compare the results to those of Benedict's 
(1989) study and note changes caused by modifications to the algorithm. In addition, 
improvements to the results derived from Mahony’s (1991) correction for atmospheric 
water absorption will be studied. Finally, the form of the Henyey-Greenstein phase 
function obtained from this study is compared to that used by Frost (1988) in order to 


make recommendations for improvement. 


Il. THEORY 


A. RADIATIVE TRANSFER 

It will be shown in the following discussion that satellite measurements of upwelling 
radiance provide a means to determine acrosol optical depth and particle size informa- 
tion. These characteristics can help determine the amount and source of aerosols. The 
satellite Wavelengths used in the approach are from the NOAA-9 AVHRR. Table 1 


gives channel bandwidths for reference. 


Table 1. AVHRR CHANNEL BANDWIDTH: (in microns) from Kidwell (1986). 
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In order to develop a practical form of the Radiative Transfer Equation (RTE) for 
use In marine aerosol calculations, it 1s necessary to review the theorv and make some 
assumptions and modifications. From Liou (1980), the form of the RTE in a plane 


parallel atmosphere for scattering radiation 1s 
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where 
L = Diffuse intensity of radiance. 
6 = Optical depth. 
w, = Single scattering albedo. 


u = cos 8 (observation zenith angle). 


it, = cos @, (solar zenith angle). 

(2 = Solid angle (8, @) (@ = azimuth angle). 

Q, = Solid angle of direct radiation singly scattered into the path. 

Q’ = Solid angle of multiply scattered radiation scattered into the path. 


p(Q, 22’) = Scattering phase function from solid angle 22 to 0’ (©). 
mf, = Incoming radiative flux (solar flux density) adjusted for variable earth-sun dis- 
tance; different for channels one and two. 

The first term on the right represents the diffuse radiance attenuated from the beam 
by single scatter and absorption while the second and third terms represent contributions 
from multiple and single scatter into the beam, respectively. 

Optical depth is defined as the amount of absorbing material lying between two 


levels and can be represented as 


H 
d -| emnaz, a 
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where o,,, 1s the extinction coefficient and H 1s the satellite height. The extinction coef- 


ext 


ficient, Which is a sum of absorption and scattering contributions, 1s defined as 
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where zr? 1s the particle cross sectional area, Q(m, /, r) 1s the extinction efficiency for 
either scattering or absorption as a function of the complex index of refraction, the 
wavelength, and the particle radius, and d\(r)/dr is the slope of the aerosol size distrib- 
ution. 

For marine particles absorption is negligible when the wavelength is Iess than one 


micron. Then o,,,, which equals the sum of the extinction from scattering and absorp- 


ext) 


tion, is nearly equal to o,.,, (Shettle and Fenn 1979). The single scattering albedo is de- 


fined as o,.,, divided by o,,, and is nearly equal to one. 


ext 
The scattering phase function, p(Q), is the distribution of scattering by particles and 
is equal to one when summed over all solid angles. More will be said about the phase 
function below. 
If multiple scattering is neglected and only single scattering of direct solar radiation 


into the path is assumed, the RTE simplifies to (Liou 1980) 


L(5; 2, 6) = L(5;; #, &) exp — (5, — 6)/H] + 
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where 6, 1s the optical depth at the bottom of the atmosphere and 6’ is the increment 
of optical depth between the bottom of the atmosphere and the level of interest. 

Additionally, over ocean areas in the red and near infrared wavebands there is very 
little radiation emitted by the earth resulting in a further simplification of the RTE at the 
top of the atmosphere to (Liou 1980) 


Oot 


L(0; u, @) = Ae 


PO) — exp[ — 0,(1/u + T/49)]}. (2.3) 
For optically thin atmospheres, where optical depth is less than about 0.1, which 1s 
true over clear ocean areas, 0,(1/u + 1/u,) < 1 and equation 2.3 reduces to (Durkee 1984) 
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L(0; 2,6) p(O)6;. (2.4) 


Now the upwelling radiance can be broken down into its various contributions as 
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where the subscripts stand for the contributions by molecules, or Rayleigh scatter, con- 
tributions from aerosols, contributions from sun glint, or specular reflection, and re- 
duction by ozone absorption, respectively. 

To get the aerosol radiance, the radiance measured by the satellite sensor was ad- 
justed by subtracting the Ravleigh scatter due to molecules following Turner (1973). 
Ravleigh scatter occurs when the size of the particle is much smaller than the wavelength 
of incident radiation and though strongly wavelength dependent, it does not vary spa- 
tially. 

An adjustment for sun glint is also necessary. Ramsev (1968) showed that sea sur- 
face reflectance is very small for red wavelengths (albedo = 0.5%) and zero for radiation 
with wavelengths greater than 0.7 microns, except where the sun-earth-satellite geometry 
results in sun glint. Therefore, geometries which could contribute to sun glint are 


avoided. Following Cox and Munk (1954), an angle of reflection is calculated from ge- 


ometry considerations. Then an elliptical area in the neighborhood of the reflection 
point where sun glint is possible is defined. This definition is based on the angular di- 
ameter of the sun. 

Since ozone absorbs weakly in channels one and two, a small ozone correction was 
added back. This correction was a modification to that of Benedict (1989). Now the 


aerosol radiance, L,, can be expressed as 
OP, 
Ly = LO; u, 6) — Ly — Lsg+ LoX ya ee (2.5) 


where 0, 1s aerosol optical depth. 

Griggs (1975) and Durkee et al. (1986) showed that there is indeed a linear re- 
lationship between aerosol radiance and optical depth. A method is now available to 
obtain aerosol optical depth from upwelling radiance by rearranging equation 2.5 to give 

Ly4u 
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By taking the ratio of channel one to channel two aerosol radiance, given by 
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a new parameter, the particle size index, is defined from which particle size information 


can be obtained: 





S\)= wy ee (2.7) 


Since for smaller particles the aerosol dependent terms, w, p(©), and 6 are larger at 
red wavelengths (Durkee 1984), smaller particles would cause S,, to be larger. The in- 
crease in optical depth can be attributed to a steeper slope in the particle size distrib- 
ution as seen from equations 2.1 and 2.2. Thus the particle size index gives an indirect 
measure of the aerosol particle size distribution. For channels one and two the radiances 


are such that S,, is normally greater than one. 


B. THE HENYEY-GREENSTEIN VARIABLE PHASE FUNCTION 
The scattering phase function describes the fraction of radiation scattered in a given 


direction. As discussed by Frost (1988), one single phase function is inadequate for 


calculating optical depth when aerosols change over global scale areas. Therefore, a 
variable phase function, the two-term Henyey-Greenstein phase function (Lenoble 


1985), was chosen and 1s given by 
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where o 1s a Weighting factor and g, and g, are asymmetry curve factors. All are empir- 


ical functions of S,, (Frost 1988) and are given by 
a= 1 2— (0.25% 55) 
g, = 1.04 — (0.08333 x S,5), 


and 


Figure 2 shows a plot of the phase function versus scattering angle for various S,, 
values. Higher values of the particle size index (smaller particles) result in a curve that 
is closer to that of Rayleigh scattering. The HG phase function uses S,, prior to cor- 
rection for Rayleigh scattering and ozone absorption. Since Rayleigh scattering, more 
in channel one, increases S,, and ozone absorption, also nore 1n channel one, decreases 
Sj. the difference is slight and the resulting particle size index and phase function do not 
vary greatly from corrected values (Durkee et al, 1991), 

The HG phase function also varies with wavelength although it is only shown here 
for channel one. This introduces an error into the calculation of optical depth for 
channel two in equation 2.6. It will be shown in the analysis of resulis how this error 


may be corrected. 


C. WATER VAPOR CORRECTION TO THE PARTICLE SIZE INDEX 

The RTE in the above calculations assumed cloud free conditions. Another impor- 
tant consideration in atmospheric absorption is the moisture present in the atmosphere. 
Water vapor can produce as much as a 15% error in the phase function and aerosol 
particle optical depth (Mahony 1991). Dalu (1986) developed an algorithm to measure 
water vapor content using the channel four and five brightness temperature difference. 
The measured water vapor content of a column of air is used to determine the “dry” 


aerosol particle size index. 
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Figure 2. The Henyey-Greenstein (HG) scattering phase function (Frost 1988) 
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In Figure 3 the effect of water vapor can be seen in the transnuttance (t) for 
AVHRR channels one and two. Note that for channel one 7 1s greater than channel two 
indicating less absorption. Channel two radiance would tend to decrease more than 
channel one leading to overestimation of S,,.. This error could be misinterpreted as a 
change in aerosol properties or it could mask real changes. Thus, accounting for water 
vapor should improve the accuracy of the results. 

Mahony (1991) used the U.S. Air Forces Low Resolution Transmittance 
(LOWTRAN)-7 atmospheric propagation model to study the effects of water vapor. 
Figure 4 shows the channel one and two relative response as Well as the resulting radi- 
ances from Rayleigh scatter, Rayleigh and aerosol scatter, and Rayleigh plus aerosol 
scatter minus water vapor absorption. The radiances increase with smaller wavelength 
as the contribution from Ravleigh scatter increases. The radiance decreases with wave- 
length in both channels but more in channel two leading to overestimation of S,,. 

An estimate of the water vapor content in an air column can be made using channels 
four and five of the AVHRR, at 1] and 12 microns respectively. Here absorption by 
water vapor is dominant while scattering 1s relatively insignificant. The benefits of using 
these channels are their collocation and temporal linkage with channels one and two, 
providing a water vapor snapshot along the same path as aerosol optical depth meas- 
urements. Although not the ideal sensor window for water vapor content, their use is 
justified within the accuracy of the approximations (Mahony 1991). 

The brightness temperature (T) is determined from the radiance received due to di- 
rect transfer from the ocean surface to the satellite sensor. Therefore the transmittance 
is directly proportional to the retrieved T. Only absorption reduces the radiance reach- 
ing the sensor, scattering is negligible. 

Dalu (1986) showed that by using the difference in brightness temperature between 


channels four and five the water vapor content could be obtained from 
w= A(T, — T,) cos 0, 


where w 1s the water vapor content in grams per square meter and T, and T, are the 
brightness temperatures in degrees Kelvin for channels four and five, respectively, and 
@ is the scanning angle. The constant A is equal to 19600 grams per degree Kelvin per 


meter squared. 
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Figure 3. Water vapor effects on NOAA-7 channels one and two (Mahony 
1991): total transmittance ts superimposed. 
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Figure 4. Water vapor influence on atmospheric radiation (Mahony 
1991): NOAA-7 response curves and LOWITRAN-7 used. 


By comparing the water vapor content to the particle size index, Mahony (1991) 
obtained the following empirical relation for the particle size index corrected for water 
vapor: 

Si 


Sn) ee 
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where w is the water vapor content in grams per square centimeter. This is the formula 


used in the calculations to obtain drv optical depth and particle size index. 
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Il. PROCEDURE 


A. THE DATA SET AND THE SATELLITE 

The in-situ. measurements were taken from the Research Vessel (R/V) 
Oceanographer, a NOAA funded research vessel, during the Radiatively Important Trace 
Species (RITS)-88 cruise conducted from 7 April to 5 May 1988. The ship’s track was 
in the central North and South Pacific Oceans from 50°N to 12°S along 170°W. Data 
measurements were obtained for every one degree of latitude and included DMS emis- 
sions in the sea and air, NSS-SO?2- aerosol concentration, aerosol size distribution, and 
total particle concentration (Charlson and Bates 1989). Optical depth was also obtained 
using a hand held photometer. 

The satellite data were collected by the AVHRR sensor onboard the NOAA-9 sat- 
ellite. See Table 1 for AVHRR channel wavelength bands. The NOAA-9 is a sun- 
synchronous, near polar-orbiting satellite with local equator crossing times of 0220, 
descending, and 1420, ascending. It has a nominal altitude of 833 km and resolution of 
1.1 x 1.1 km at sub-point. 

The satellite data utilized were for the same time frame and area as the ship’s track 
from the digital archives of the National Environmental Satellite and Information Ser- 
vice (NESDIS). The 57 passes were stored as Global Area Coverage (GAC) data; GAC 
is an archiving technique that also reduces the resolution at sub-point to about 4.0 x 4.0 


kim. 


B. PROCESSING FACILITIES 

The data were processed at the Interactive Digital Environmental (IDEA) Labora- 
tory at the Naval Postgraduate School in Monterey, California. The main computer 
was a Digital Equipment Corporation VAX 8250 with graphics display terminals. 
Processing included loading and reading the magnetic tapes to memory and then proc- 
essing the data. The multi-channel analysis routine used to analyze the data was a 
FORTRAN program run on the VAX computer. Visual images were also produced in 


the IDEA lab and helped in the interpretation and analysis of the data. 


C. THE ALGORITHM 
The satellite data were analyzed using a scheme developed by Pfeil (1986) and 


modified by Frost (1988). In addition, corrections were made to solar irradiance to ac- 





count for variations in earth-sun distance, for over correction of ozone absorption, and 
to correct the computer code for a mathematical error. Calculations were also made 
using the water vapor correction to S,, described by Mahony (1991). 

The pixels were subjected to a number of filters, described below, and composited 
in one degree by one degree boxes. To study the relationship between satellite and 
shipboard measurements, channel one aerosol optical depth and the particle size index 
were further filtered geographically to within five degrees longitude of 170° W and plus 
or minus three days of the ship’s track. 

The first step in the algorithm is to test for valid data. After reading in the data and 
determining their geographical position, only data equatorward of 70°N or S is accepted. 
In addition, the data are required to be over ocean areas; here upwelling radiance from 
the surface 1s minimized. Next, sun angle is checked. Solar zenith angle must be less 
than 70° to avoid longer and more complicated optical paths. Data with specific ge- 
ometries 1n relation to the satellite are also eliminated to avoid sun glint (Cox and Munk 
1954). This eliminated a large portion of the western half of each pass. 

Before testing the pixels further, an adjustment to channels one and two were made 
to account for ozone absorption. This is of greatest concern for a wavelength of 0.69 
mrenons (Durkee et al. 1991). 

The calibrated pixels are then further processed to determine cloud contamination. 
If the brightness temperature of channel four is less than 273°K, the pixel is assumed to 
be contaminated by high clouds. Next, if the channel two albedo 1s greater than 40%, 
the pixel 1s assumed to be contaminated bv low clouds. Then, if the channel one divided 
by channel two albedo ratio 1s greater than 1.5, the pixel 1s considered clear since cloud 
scattering 1s weakly dependent on wavelength. If the ratio is less than 1.5, the pixel 
undergoes further testing. The channel two albedo 1s compared to that of the pixels 
above, below, and to the left and right. If the difference between the maximum and 
minimum counts ts less than five, the pixel is clear, otherwise it 1s partly cloudy. A count 
iS approximately 0.1069% albedo (Kidwell 1986). 

Once a clear pixel has been determined, the Rayleigh radiance 1s subtracted from 
channels one and two (Turner 1973) and the optical depth and particle size index are 


calculated using equations 2.6 and 2.7. The derived parameters are: 


The Spatial Variables. Channel one and two radiances were used to obtain optical 
depth and particle size index. The “dry” values were also obtained using the algorithm 
of Mahony (1991). This information was composited temporally over the entire data 
set and was used to create images. 
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The Geometry Variables. These are variables such as normalized measured radiance 
as a function of scattering angle, latitude, and particle size index. Thev were com- 
posited geographically across all longitudes as well as temporally. One parameter, the 
normalized measured radiance as a function of scattering angle for different particle 
size indices, was used to study the two-term HG phase function. 


The Daily Variables. The particle size index and channel one optical depth both un- 
corrected and corrected for water vapor filtered to be within five degrees of longitude 
and plus or minus three days of the ship’s track. This information was used for com- 
parisons with in-situ measurements. 


The calculated values were then stored in data files and used as needed to produce 
the necessary images and plots for comparison with those of Benedict (1989) and the 


measurements of the R/V Oceanographer. 


IV. RESULTS AND DISCUSSION 


The results are organized to show improvement to the analysis done by Benedict 
(1989), to show improvements to aerosol characteristics based on corrections for column 
water vapor as parameterized by Mahony (1991), and to investigate improvements to the 
HG variable phase function. 

Before looking at the results, it is useful to examine a summary of the synoptic sit- 
uation during the 7 April - 5 May 1988 time frame. The synoptic situation is important 
to the discussion since it helps identify sources of aerosol particles. Knowing the sources 
of the aerosols helps in interpreting the results. 

Representative Navy Operational Global Atmospheric Prediction System 
(NOGAPS) surface pressure and 500mb height patterns are shown in Figures 5 and 6. 
A generally westerly flow exists over the Asian continent tn the 35-45°N latitude band 
indicating a source for natural and anthropogenic continental aerosols. Natural aerosols 
would tend to be larger in size than the marine aerosols being studied. Benedict (1989) 
discusses the occurrence of KOSA, or dust clouds from storms in the Gobi desert during 
April. Given typical translational rates (Shaw 1980), the dust from these storms could 
influence the aerosol concentration in the central Pacific. 

Another environmental factor discussed by Benedict (1989) was the eruption of the 
Kilauea volcano in Hawai (17°N, 157°W) producing a continual source of aerosols 
during this time frame. These aerosols Were initially dominated by gaseous sulfur com- 
pounds which later changed to the solid phase in the atmosphere and grew as they mi- 
grated downstream. The debris from the eruption was carried across the path of the R/V 
Oceanographer’s track by the Northeast Tradewinds and influenced the aerosol meas- 
urements. 

The composite optical depths for the entire experimental period are shown in Fig- 
ures 7 and §. The optical depth is an indication of both the number and size of aerosols. 
The uncorrected optical depth is presented in Figure 7 which shows an area of relatively 
low optical depth in its lower right corner. A higher optical depth was expected due to 
the volcanic aerosols. The reason for the low optical depth is because most of the 
volcanic enussions are gaseous and very small and therefore do not greatly effect optical 
depth. It is also possible that the parameterization of the phase function is not correct 


leading to errors in optical depth. Figure 7 shows an area of high optical depth at the 
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NOGAPS 1200 UTC 10 April 1988 Surface Analysis (Benedict 1989) 


Figure 5. 





NOGAPS 1200 UTC 10 April 1988 500mb Analysis (Benedict 1989) 


Figure 6. 
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western edge of the area which extends eastward across the upper third of the figure. 
This indicates a potentially large number of particulates which could be either 
anthropogenic or natural in origin. When the particle size index is examined below, 
further information about the source of the aerosols may be inferred. From Figure 2 it 
is seen that for a given radiance, when the particle size index decreases the phase func- 
tion decreases. This leads to an increase in optical depth for most scattering angles. This 
is the case for the corrected optical depth in Figure 8. Here again a high optical depth 
1s noted on the left periphery of the figure. It extends eastward into the upper and cen- 
tral areas of the figure. This mav be due to dust intrusions from Asia. Another possi- 
bility for the increasing optical depth to the west 1s related to the retrieval process. This 
is the result of the satellite sensor’s aperture moving to the left and to the right of the 
of the satellite's path to view side regions. The satellite is looking mainly at backscatter 
When it looks directly at the Earth. As the sensor looks to the side it 1s no longer viewing 
backscatter but smaller scattering angles. This results in a decrease 1n the phase function 
and a corresponding increase in the optical depth. This is an area that should be inves- 
tigated further in future work. 

Figure 9 shows Benedict’s (1989) composite optical depth averaged geographically 
for all longitudes. There is a peak in the vicinity of 40°N indicating a potentially greater 
number of particles in the region. This could be a reflection of the advection of Asian 
particulates into the area. These could be natural dust or anthropogenic pollutants. 
There 1s another peak near 8°N. Since Figure § revealed little optical depth from the 
erupuion of Kilauea, this peak is probably due to DMS production in the tropics. 

The optical depths obtained from this study, shown together in Figure 10, reveal a 
profile with greater structure, especially in lower latitudes where there is a peak at 
22°N not revealed by Benedict (1989). The peak at 6°N is also shifted north to near 10 
°N. This shows possible error in the HG phase function parameterization. The increase 
in optical depth south of 20°S is subject to error as it is near the data set boundary. The 
values of optical depth are larger than those of Benedict due to the ozone correction. 
This correction increases channel one radiance resulting in a higher optical depth. 

Figures 11 and 12 present the composite particle size index and particle size index 
corrected for water vapor for the entire experimental period. The uncorrected particle 
size index shows a plume in red extending westward from the Kilauea volcano. The 
higher particle size index indicates smaller particles. To the north in the midlatitudes 1s 
a darker area of blue indicating larger particles. This could be attributed to dust from 


the Asian continent being advected to the east. The large area of dark blue to the 
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Figure 7. Composite channel one optical depth: the area extends from roughl 
"N to 12°S and 162°E to 142° W 
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Figure 8 Dry composite channel one optical depth: — the area extends from roughly 
BG tomes exci 162 to 142° W. 
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Figure 9. Composite channel one optical depth (Benedict 1989) 
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Figure 10. Composite uncorrected and dry channel one optical depth: the uncor- 
rected values are smaller than the dry values. 
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western side of the image is unidentified although it could represent more dust. The 
Somewhat lower particle size index on the left side of the figure might also be related to 
the retrieval process discussed above. 

The particle size index corrected for water vapor, Figure 12, shows lower values. 
Since water vapor causes an overestimation of S,,, this was expected. Now the plume 
from Kilauea is less evident but the dust influence to the north is more significant. 

Figure 13 shows the composite particle size index obtained by Benedict (1989). The 
results are averaged geographically for all longitudes. A peak in the neighborhood of 
30-40°N is evident indicating possible influence by smaller particles being advected from 
the west. Another peak is noticed at about 6°N indicating smaller particles, possibly a 
result of DMS production in the tropics. 

The plots of composite particle size index and dry particle size index for this study 
are both shown 1n Figure 14 for cemparison. They show more structure and definition 
than the plot of Benedict (1989). The values are higher than those of Benedict (1989) 
due to the modifications in the algorithm. In Benedict’s (1989) study ozone was over 
corrected for leading to a proportionately higher radiance in channel two and thus lower 
particle size indices. The uncorrected index is higher than the corrected index as before. 
Both curves show the nudlatitude peak of Benedict (1989) but also show more of a peak 
in the tropical regions. The Benedict (1989) plot decreases toward the equator while the 
plots of this study increase. This supports the changes to the algorithm since the particle 
size index is expected to increase in the tropics because of higher DMS production. In 
addition, the dry index has a larger correction for the tropical regions due to the greater 
effect of water vapor there. The increase in S,, near 60°N may be due to a reduction of 
valid data near the edge of the area. 

The next data sets are in an area of plus or minus five degrees of longitude from the 
ship’s track and within plus or minus three days of the ship’s position which limits the 
data used in the north-south direction. This is to get a better comparison with the 1n- 
situ measurements. 

Optical depth near the ship’s track from Benedict (1989) is shown in Figure 15. A 
large double peak at about 40°N is evident with a smaller peak at about 33°N. In Figure 
16 this study shows the same peaks at midlatitudes but also a large peak at 22°N. This 
is most likely due to volcanic particles which may have been averaged out in the previous 
composite (Figure 10) which included the entire data set. The corrected optical depths 
on the same plot are even larger. Higher optical depths were expected due to both the 


Water vapor correction to the particle size index and the ozone correction. The optical 
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Figure 13. Composite particle size index (Benedict 1989) 
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Figure 14. Composite uncorrected and dry particle size index: The uncorrected 
values are larger than the corrected values. 4 
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depth measured by a hand-held sun photometer on the ship 1s also plotted on Figures 
15 and 16. Note that optical depth is scaled differently in the two figures. The magni- 
tude of the optical depths obtained in this study appear to be supported by the 
photometer measurements. This indicates support for the corrections made to the al- 
gorithm. There appears to be a northward shift of about 5-7° in the satellite estimated 
optical depth compared to the ship measurements. This could be caused in part by the 
wav the data 1s averaged by three days around the ship’s track. It 1s possible that sat- 
ellite data within one day of the ship’s track would produce better results. The reason 
this wasn’t done was to ensure enough data were available to fill in times when cloudy 
pixels were thrown out. A second reason for the shift may be due to inaccuracies in the 
positioning of either the ship or satellite data and a third alternative is that the 
parameterization of the HG phase function 1s incorrect. Changes to improve the phase 
function are discussed below. 

In-situ measurements of condensation nuclei are shown in Figure 17 although they 
too are shifted several degrees south of the satellite data. The condensation nuclei are 
not the same as CCN, CCN are a subset of the condensation nuclei, being the 
condensation nuclei capable of nucleating with water vapor to form cloud droplets. 
However, for the purpose of optical depth measurement they can be treated similarly. 
Another peak not resolved by Benedict’s (1989) or this data set is evident at 6 and 
ase 

Figure 18 from Benedict (1989) shows the particle size index in the neighborhood 
of the ship’s position. It is a little more definitive than the previous composite and 
reveals a midlatitude peak around 30-40°N, another at about 20°N, and another near 
10°N. A distinction can be made between an influence from Kilauea volcano at 20°N 
and DMS production leading to non-sea-salt sulfate aerosols near the equator. 

Figure 19 shows the combined uncorrected and drv particle size index near the re- 
search vessel's track for this study. The peaks are more pronounced than Benedict 
(1989) and there is a minimum at 22°N. where Benedict (1989) shows a peak. The values 
of S,, are higher overall and the difference between uncorrected and dry particle size in- 
dex is larger near the equator. There is also an increase in the results moving south of 
the equator as well as a peak not seen in the Benedict (1989) results, possibly due to 
DMS production. 

Next, the non-sea-salt sulfate data obtained from the R, V Oceanographer are plotted 
in Figure 20. Tlere there is a peak clearly evident centered at about 20°N. This is strong 


indication of smaller particles, most likely from the volcanic fallout. These results point 
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Figure 15. Channel one optical depth along ship’s track (Benedict 
1989): photometer measurements superimposed. 
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Figure 16. Uncorrected and dry channel one optical depth along = ship’s 
track: photometer measurements superimposed 
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Figure 17. Condensation nuclei along ship’s track (Benedict 1989) 
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Figure 18. Particle size index along ship’s track (Benedict 1989) 
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Figure 19. Uncorrected and dry particle size index along ship’s track: the uncor- 
rected values are larger than the dry values. 


37 


LATITUDE PROFILE 


60°N 
SO°N 
40°N 
SO? N 
20° N 
1C°N 
EQ 

1s 


20°S 
0 S 10 ts <0 ao 30 J5 40 


NSS-SO, (nm-m~7>) 


Figure 20. Non-sea-salt sulfate along ship’s track (Benedict 1°8°) 


to errors in the phase function parameterization since there 1s no corresponding peak in 
particle size index as would be expected. 

The next plot to be studied is the geometry variable used to examine the correlation 
between optical depth and phase function. By plotting the normalized measured radi- 
ance as a function of scattering angle for various particle size indices, any correlation 
between the phase function and particle size index is separated out. This can be seen 
by rearranging equation 2.4 to give 

ee ais. 


Wok, 





The left hand side represents values that are all measured or known while the two terms 
on the right must be parameterized (the phase function) or calculated (the optical 
depth). By using the optical depth that best fits the data for a given measured particle 
size index, the phase function can be adjusted and improved. The phase function used 
by Frost (1988) and Benedict (1989) was simply one way, not necessaruy the best way 
for this data, of parameterizing the phase function. 

The results of P versus © are shown in Figure 21 and can be compared to Figure 2, 
the FIG phase function. To make an accurate comparison, the IIG phase function used 
bv Frost (1988) should be multiplled by a mean optical depth. This is done for two 
values of the particle size index in Figure 22 and allows for a direct comparison between 
the HG phase function and study results. For a particle size index of 1.7, the 11G phase 
function was multiplied by a mean optical depth of 0.2 and for a particle size index of 
2.1, the HG phase function was multiplied by a mean optical depth of 0.1. These values 
of mean optical depth are only approximate. Further evaluation of the study results 
would be required for more accurate values. 

The shapes of both curves are similar but in the area of 0-20° scattering angle the 
results of this study are up to three orders of magnitude lower. This is due partly to the 
observations 1n the 0-60° bin. Most of the pixels in this bin have a scattering angle be- 
tween 50 and 60° resulting in a P weighted toward the lower values of the higher scat- 
tering angles. The total number of pixels evaluated for this bin was also less than the 
other bins potentially contributing to an unrepresentative P in this area. For the particle 
size of 1.7 in this study, there also appears to be a shift in minimum P values toward 


higher scattering angles than for the H1G phase function. 


39 


By adjusting phase function parameters a better fit can be obtained. In Figure 23, 
after multiplying the HG phase function by a mean optical depth of 0.15 and using the 
Si. = 1.7 data curve, the following values from equation 2.8 were obtained: a was set 
to 0.95, g, and g, were set to 0.5. In this way the form of a, g,, and g, can be determined. 

Another aspect of the problem 1s that the results were always obtained using the 
HG phase function parameterization for channel one. Channel two results are also 
available but they are computed using the phase function for channel one. The phase 
function is wavelength dependent and could result in some error for channel two calcu- 
lations. Using the same procedure as above but for channel two, it should be possible 
to derive a channel two correction to the Henyey-Greenstein phase function. This would 
result in knowing all parameters in the simplified RTE except optical depth, which could 
then be calculated directly for either channel. A new particle size ratio of optical depths 


from channels one and two would then more directly measure particle size information. 
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Normalized radiance (P) versus scattering angle (Q): for various S,,. 
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Scaled HG phase function and P versus ©: for particle size index of 
1.7 (solid line) and 2.1. 
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Figure 23. Scaled and adjusted HG phase function and P versus ©: for a particle 
size index of 1.7. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This study examined an algorithm that extracts aerosol particle size index and op- 
tical depth from satellite observations of upwelling radiance over clear ocean areas. By 
utilizing the difference in sensitivity between channels, particle size information is 
measured and optical depth can then be calculated. Data retrieved from the NOAA-9 
AVHRR sensor was reanalyzed and compared to shipboard measurements made by the 
R/V Oceanographer during the RITS-88 cruise. In this fashion the satellite retrieval 
technique 1s verified and our understanding of marine aerosols on a regional scale is in- 
creased. 

This study improved upon the results of Benedict (1989) after making corrections to 
account for variable earth-sun distance on solar irradiance, making an adjustment for a 
previous over correction to ozone absorption and correcting the computer code for a 
mathematical error. These combined corrections had a greater effect on the results than 
anticipated. The direction of the corrections was positive, however. An improvement 
to the algorithm, which takes into account the absorption of column water vapor, was 
examined and also found to improve the results. 

Finally, the correlation between aerosol optical depth and the Henyey-Greenstein 
phase function was examined. It was concluded from the plots of P versus © for various 
particle size indices that some correlation did exist. An adjustment to one case of the 
HG phase function was made with excellent results. It 1s recommended to calculate a 
new set of phase functions that fit the experimental data. This would account for the 
correlation between optical depth and scattering angle thus making optical depth and 
phase function independent of one another. Then optical depth could easily be calcu- 
lated and a ratio of optical depths would be a new, more accurate particle size index. 
This 1s because the slope of the particle size distribution, which 1s part of optical depth, 
would be more directly compared. 

Another recommendation of this study is to separatelv calculate a correction to the 
HG phase function for channel two. Then the channel two phase function could be used 
in the calculations of channel two optical depth instead of the channel one phase func- 
tion. It is also recommended that the simplified RTE be adjusted for multiple scatter. 
The single scattering approximation is good up to an optical depth of about 0.1. The 


results obtained show that the optical depths were higher than 0.1 leading to errors in 
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the results. Other regional studies should be undertaken to continue examining and 
improving the algorithm and its performance. Global scale studies are also useful but 
it is less likely that there will be in-situ data available for verification. 

While this procedure appeared effective in revealing variations in aerosols in the re- 
gion of interest for known sources, it appeared to have only moderate sensituvity to the 
NSS-SO32-’s of interest. This makes studying DMS and its role in global chmate regu- 
Jation more difficult. It is therefore important that improvements to the algorithm 
continue to be made. As discussed by Benedict (1989), studving areas where there is less 
interference from natural or human-made sources is one possibility as is studying regions 
of higher phytoplankton activity, such as the tropics. Additional satelhte-inferred aer- 
osol studies are required to guide future algorithm improvements and to provide crit- 


ically needed measurements of global aerosol behavior. 
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